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Abstract 

The genus g free energies of matrix models can be promoted to modular invariant, non- 
holomorphic amplitudes which only depend on the geometry of the classical spectral curve. 
We show that these non-holomorphic amplitudes satisfy the holomorphic anomaly equations 
of Bershadsky, Cecotti, Ooguri and Vafa. We derive as well holomorphic anomaly equations 
for the open string sector. These results provide evidence at all genera for the Dijkgraaf-Vafa 
conjecture relating matrix models to type B topological strings on certain local Calabi-Yau 
threefolds. 
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1 Introduction and conclusions 

Topological string theory has been a fascinating laboratory to explore issues in string the- 
ory with important connections to other branches of physics and mathematics. The basic 
problem in topological string theory is to compute and understand closed and open string 
amplitudes on different geometric backgrounds. Particularly important among these are 
Calabi-Yau (CY) manifolds. Different techniques have been developed for the computation 
of these amplitudes. For type B topological strings on CY manifolds, a powerful method to 
solve the closed sector of the model are the holomorphic anomaly equations of [6]. These 
equations control the t-dependence of the closed string amplitudes F^^\t,t), and when com- 
bined with extra boundary conditions, they lead to explicit answers, see [T9l [T7] for recent 
progress in this direction. 

On the other hand, in some special backgrounds one can use large N dualities and 
geometric transitions to compute the holomorphic limit of F^^\t,t), which will be denoted 
by F^^\t). In a groundbreaking paper [11], Dijkgraaf and Vafa conjectured that on certain 
noncompact Calabi-Yau manifolds, where the geometry reduces to a complex curve, the 
F^^\t) are given by the genus g free energies of a matrix model in the 1/A^ expansion 
(see pl| for a review). The curve encoding the CY geometry is then identified as the 
classical spectral curve of the matrix model. Given the connection between topological strings 
and type II superstrings, this has made possible to compute certain protected quantities in 
super symmetric gauge theories by using matrix model technology. The connection between 
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topological strings on noncompact CY manifolds and matrix models has been extended 
recently to the mirrors of toric geometries |25j . 

The conjecture of Dijkgraaf and Vafa was verified at the planar level in [11], and at genus 
one in [22llT2j. In the simple case of the cubic matrix model, evidence for the conjecture was 
given at genus two in [1^, where it was shown that the matrix model expression for 
is the holomorphic limit of a particular solution to the holomorphic anomaly equations for 
the corresponding local curve. 

In this paper we will give further evidence for the Dijkgraaf- Vafa conjecture and the 
related results of [25j, by using recent advances in the solution of matrix models at all orders 
in the 1/A^ expansion obtained in [131 [151 [9] ciiid more recently in \X6l- One of the outcomes of 
these advances is that, as explained in [16], one can extrapolate the matrix model procedure 
and define a series of amplitudes F^^\t) and correlation functions Wg{pk) for any algebraic 
curve H{x, y) = 0. When this curve is the classical spectral curve of a matrix model, one 
obtains in this way the 1/A^ expansion of the free energies and correlation functions. Using 
the results of [16j, we will show that it is possible to construct non-holomorphic free energies 
F^^\t,t) and correlation functions which satisfy the following conditions: 

• In the holomorphic limit t — > oo they reduce to the holomorphic amplitudes F^^\t) 
associated to the spectral curve. 

• They are invariant with respect to the symplectic modular group of the curve. This is 
in contrast to their holomorphic limit, which does not have good modular properties 

• They satisfy the holomorphic anomaly equations of [6]. 

This gives a procedure to obtain the full non-holomorphic couplings F^^\t, t) in the con- 
text of matrix models, by using the requirement of modular invariance as in [T], [T0|, and 
shows that these couplings obey the equations of [B] that characterize topological string am- 
plitudes. This implies in particular that the matrix model free energy at genus g (extended 
non-holomorphically in this way) must be equal to the type B free energy of the noncompact 
Calabi-Yau manifolds considered in [HI [25], up to a holomorphic modular invariant quan- 
tity. Proving the Dijkgraaf- Vafa conjecture reduces now to proving that these holomorphic 
invariant functions vanish at all genera. 

An interesting spinoff of our work is that we can determine the holomorphic anomaly 
equations for open string amplitudes with respect to the closed moduli i, at least in the 
local case. A direct string theory derivation of these equations has shown to be rather 
elusive [H H] , and we expect them to be useful in future investigations of the open sector in 
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the general case. The equations we find turn out to be closely related to the holomorphic 
anomaly equation for correlation functions obtained in [6], as we will explain in section 4. 

This work has other interesting implications. If one considers the hyperelliptic curves 
that appear in Seiberg-Witten theory as classical spectral curves of a matrix model, the 
amplitudes F^^^ constructed with the recursive procedure of [16] should be closely related 
to the gravitational couplings of Seiberg-Witten theory introduced by Nekrasov in [27J. It 
is known that these gravitational couplings can be promoted to non-holomorphic quantities 
that satisfy the holomorphic anomaly equations associated to the Seiberg-Witten curve 
[m [IT]. It follows from the results in this paper that these gravitational couplings differ 
from the F^^'^ of the spectral curve at most in a holomorphic, modular invariant object, and 
it is natural to conjecture that they are in fact equal. The existence of a matrix model 
description of Nekrasov's Fg has been proposed in [7]. 

Our work also leaves some open problems. First of all, although it goes a long way 
towards proving the conjecture in [TT], it would be nice to finish the proof by fixing the 
remaining ambiguity. At a more conceptual level, the t dependence of the amplitudes has 
been introduced, following [16], by requiring modular invariance and modifying the recursive 
construction of [16] accordingly. It would be very interesting to find a more intrinsic way to 
introduce this non-holomorphic dependence in the theories which are defined by matrix inte- 
grals. This might shed some light on the problem of background independence in topological 
string theory and of the wavefunction behavior of the topological string partition function 

m- 

The organization of this paper is as follows. In section 2 we review the holomorphic 
anomaly equations and the connections between topological strings and matrix models. Most 
of the technology needed in this paper has been developed in [16j, and section 3 contains 
some results of that paper which will be essential in our derivation. Section 4 contains the 
derivation of the holomorphic anomaly equations. We actually provide a direct derivation 
and also a combinatorial derivation which makes contact with the recent results of [1]. 

2 Holomorphic anomaly and topological strings 
2.1 The holomorphic anomaly equations 

Topological string theory on Calabi-Yau threefolds (see [^ [2U] ) for recent reviews) is 
defined by coupling the twisted J\f = 2 sigma model to topological gravity and it comes in 
two versions: the A model (related to Gromov-Witten theory) and the B-model (related to 
deformations of complex structures). The topological amplitudes depend on a set of moduli 
that can be regarded as marginal deformations of the underlying M = 2 theory. If we denote 
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the chiral and antichiral operators of the theory that correspond to marginal directions by 
0/, 0j, J = 1, ■ ■ ■ , n, we can deform the action in the standard way by using their descendants 



S^=2 = S{0) + J2t^ # + / (2.1) 



1 = 1 •'^9 1 = 1 

In the closed string sector we are interested in computing the genus g free energy F^^^ as 
well as the correlation functions of integrated chiral operators, 

Ct-..={l 4!S-f 4f?l) (2.2) 

Since the t perturbation of the action can be written as a BRST-exact term with respect to 
the topological symmetry of the theory, one expects F^^^ and the correlation functions to be 
i-independent. However, as discovered in [5l [6], this is not the case once the twisted J\f = 2 
sigma model is coupled to gravity. The BRST-exact terms give a boundary contribution 
that can be evaluated recursively, and one finds holomorphic anomaly equations for the t 
dependence of the amphtudes. 

In order to write down the holomorphic anomaly equations, we need some basic ingre- 
dients from special geometry. The Zamolodchikov metric on the moduli space of marginal 
deformations is Kahler and can be derived from a Kahler function 

Gjj = djdjK. (2.3) 

The amplitude at genus zero Fq, also called the prepotential, leads to a holomorphic three 
point function 

The holomorphic two-point coupling is given by 

and plays an important role, as we will see later. Out of these quantities one can define a 
tensor _ _ 

C^' = e^^G^^G^^^TO- (2.6) 

The Zamolodchikov metric (12.31) defines a connection on moduli space with Christoffel sym- 
bols _ _ _ 

Tfj = G^''d,Gjj„ Tfj = G''^djG,a. (2.7) 
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It turns out that the F^^) 

are sections of a hne bundle on the moduh space, and the 

covariant derivative acting on them is given by 

Di = di-Ti + {2- 2g)djK, (2.8) 

where K is the Kahler function appearing in (I2.3p . With these ingredients, the holomorphic 
anomaly equations of [6] can be written as 

^F(f) = - CLJ {DiDjF,.^ + DiFhDjFg.f,^ , g>2. (2.9) 

h=l 

The equation for g = 1 is slightly different and reads 



didjF^ = \CjklC^^ - - l]Grj. (2.10) 



Of course, the holomorphic anomaly equations determine F*^^^ (t, t) only up to a holomorphic 
function of the moduli. This integration constant is called the holomorphic ambiguity. The 
holomorphic limit of F'^^\t,t) is defined as 

F(5)(t) =_lim F(^)(t,t). (2.11) 

This is the limit which is for example appropriate to make contact with Gromov-Witten 
theory, in the A model. It has been observed [lO, H] that the non-holomorphic couplings 
F^^\t,t) are invariant under the symmetry group acting naturally on the moduli space. In 
contrast, the holomorphic limit (12. lip breaks the modularity properties. 

We finally point out that the correlation functions (12.21) can be written as covariant 
derivatives of the F^^\t,t): 

C\f...,^ = Dj,...Di^F^^\ (2.12) 
and they satisfy the holomorphic anomaly equation 

9 n 



n 

-(2^-2 + n-l)5^G^,Ci^..,_,^^^...,„. 



s=l 
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2.2 The local Calabi— Yau case and the Dijkgraaf— Vafa conjecture 

A particularly interesting class of topological strings are those defined on non-compact 
Calabi- Yau manifolds. These are commonly referred to as local Calabi- Yau manifolds. 
In the local B model the Calabi- Yau geometry is encoded in a Riemann surface S„ of genus 
n [21] (this should not be confused with the worldsheet of the string), and the resulting 
structure is significantly simpler than in the global case. The number of moduli is given by 
n, the genus of the curve, and they can be parametrized by 

cc / y{x)dx, I = l,---,n, (2.14) 

J A' 

where are the compact A-cycles of S^, and y{x)dx is a meromorphic one- form on S„ 
which can be constructed from the equation for the spectral curve. The proportionality 
constant in (12.141) will be determined below by a detailed comparison with the matrix model 
results. The imaginary part of the two-point coupling (12.51) is now a positive definite matrix, 
and the Zamolodchikov metric is simply given by 

Gjj = -i{T-f)jj. (2.15) 

As noticed in [231 [HI [I]! the holomorphic anomaly equations also simplify in the local case. 
First of all, the tensor entering in (12.91) is now simply given by 

C^' = G^^G'^Cjnm = -[(r - frY'^r - r)-^^^™. (2.16) 

The covariant derivative acting on the F^^^ amplitudes involves now only the Christoffel 
symbol (i.e. the term involving djK is no longer present). Notice that the Christoffel 
symbol can be written as 

Tfj=[iT-f)-r'^djTjM, (2.17) 

With these modifications, the holomorphic anomaly for F^^\ g > 2 is still given by (12. 9p . 
For g = 1 the last term of (I2.10p is not present. The holomorphic anomaly equations for the 
correlation functions are given by (I2.13P but the contribution of the last line is absent in the 
local case. 

Since it will be useful later one, we now write the holomorphic anomaly equations for 
g >2, in the local case, in terms of a generating functional. We introduce. 



oo 

exp^ 

9=0 



oo 

^Y^N^-^3F^3\t)j, (2.18) 
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where plays the role of a genus-counting parameter. We also define, 

Z = exp{-A^2^(°) - Z = exp{-A^^F(°)}Z. (2.19) 

For g > 2 we want to express the antiholomorphic derivative of Z in terms of Z. Using that 

dl log Z = jdljZ - ^diZdjZ (2.20) 

we can write (12.91) in the local 

ldj,Z = -^lc^'[dhZ -Tfjdj,z). (2.21) 

In the local case, when the Calabi-Yau geometry reduces to a Riemann surface S„ of 
genus n, the natural group acting on the moduli space is the mapping class group, and there 
is an induced action of the symplectic modular group 

Sp(2n,Z) (2.22) 

on the topological string amplitudes F^^\t,t) and their holomorphic limit. The non- holo- 
morphic amplitudes F^^\t,t) of the local geometry turn out to be modular invariant under 
the action of this group [HI [H [lO], while the F^^\t) do not have good modular properties. 

An interesting class of local Calabi-Yau backgrounds are the ones described by the equa- 
tion 

uv = H{x,y), H{x,y)=y^-iW'ix)f + fix) (2.23) 

where W{x), f{x) are polynomials of degree d + 1, d — 1 respectively. The Riemann surface 
H{x, y) = associated to this geometry is the hyperelliptic curve 

y' = {W'{x)f - fix) (2.24) 

of genus n = d — 1. Dijkgraaf and Vafa conjectured in [TT] (see [21] for a review) that the 
holomorphic amplitudes F^^\t) of type B topological string theory on the backgrounds of 
the form (12.231) are given by the 1/A^ expansion of a matrix model with spectral curve (I2.24p . 
This is in fact a matrix model with potential W{x) and generically with d — 1 cuts around 
its extrema. This conjecture has been verified at the planar level in [11], at genus one in 
[221 [12], and at genus two in [IB]. 

Another interesting class of local geometries is given by the mirrors of toric Calabi-Yau 
manifolds. In this case, the Riemann surface is a complex curve in "exponentiated" variables 
which can be written in the form 

I 

y^ = M\x)l[{x-x,), (2.25) 

1=1 
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where M{x) is in this case a transcendental function. As shown in [16] and reviewed below, 
one can generalize the matrix model idea and associate a series of free energies F^^\t) 
and correlation functions to any algebraic curve. It was recently proposed in [25] that 
holomorphic open and closed type B topological string amplitudes for the mirrors of toric 
background are given by the correlation functions and free energies, respectively, associated 
to the spectral curve fl2.25p with the procedure of \i3\. 

We then see that, for local CY backgrounds, where the geometry reduces to a complex 
curve, the holomorphic limit of topological string amplitudes in the B model has been con- 
jectured to be given by the free energies and correlation functions associated to the complex 
curve by the recursive procedure of [16]. In the example (12.231) this coincides with the 1/A^ 
expansion of the free energy and correlators of a matrix model with potential W{x). In the 
rest of this paper we will give strong support to this conjecture by showing that the quan- 
tities defined in [16] can be extended to non-holomorphic objects which have good modular 
properties and satisfy the holomorphic anomaly equation fl2.9p in the local case. 



3 Review of matrix models 

3.1 Formal matrix models and algebraic geometry 

Random matrix models became very popular in the 80's and 90's after it was discovered by 
't Hooft [28] and then in [8] that the Feynman graphs of Hermitian matrix integrals are in 
fact discrete surfaces. Formal random matrix models are thus generating functions which 
count discrete maps of given topology. For instance the 1 matrix integral: 

Z,MM = y'rfMe-fT'^^W (3.1) 

counts discrete maps made of triangles, squares, ... fc-gones where k < deg\^. The 2 matrix 
integral 

Z2MM = j dM,dM2 e-f Tr(yi(MO+y.(M,)-M,M,) (3 2) 

counts discrete maps made of polygons of two possible colors (or call it two possible spins 
±). It was introduced by Kazakov [20j as an Ising model on a random map. 

In all cases, the formal integral is defined as a formal series in its small t expansion, and 
it turns out that to any order in t, one has a "topological expansion": 

oo 

Fmm = In ZuM = N'-'^F^^^ (3.3) 

5=0 
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where F^^^ is the generating function which counts maps of genus g. (see [H] for a review 
on this subject). 

The goal was then to compute exphcitly the coefficients F^^^ in that expansion. Many 
methods have been invented, but the only one which has really succeeded so far is the 
so-called loop equation method which just amounts to integrate by parts or equivalently 
write Ward identities or Virasoro (or W-algebra) constraints. It was soon realized that 
loop equations allow in principle to find F^^^ for any g [3, 2j, but a more explicit solution 
(which can be represented diagrammatically) was found only recently first for the 1-matrix 
model [T3], then for the 2-matrix model [l5l[9], and more recently for the matrix model with 
external field [16j. 

Then it was understood in [16] , that the solution of matrix models loop equation extends 
beyond the context of matrix models and is in fact a property of algebraic geometry in 
general. The diagrammatic method of [IE] allows to compute explicitly all the free energies 
F^^^ appearing in (13. 3p as well as all correlation functions (generating functions of discrete 
open surfaces of genus g with k boundaries): 

where the mean value refers to the measure in (13.11) and the variables Xj are fugacities 
associated to the lengths of boundaries. 

Let us now give a brief summary of the construction of |16j . 

Given an arbitrary algebraic curve H{x, y) = 0, of genus n, one constructs recursively a 
sequence of multilinear symmetric meromorphic forms wjf^ and some scalars F^^'^'s in terms 
of residues localized at branchpoints of the curve. 

Many properties of those wjf^^s and F^^^^s are described in [16], in particular variations 
of conformal structure, integrable structure, invariance under symplectic transformations of 
the curve, homogeneity, and the one which interests us now is modular transformations. 

Indeed, the construction of [16] is based on the Bergmann kernel, which is defined for a 
given canonical basis of cycles. When one changes the basis of cycles, the Bergmann kernel 
changes, and the wj^^^^s and F^^^^s change accordingly. In [T6], a "modified" Bergmann 
kernel was introduced in order to easily take into account those modular changes. The 
modified Bergmann kernel depends on an arbitrary complex symmetric matrix k, and a 
modular change of cycles merely amounts (as far as only the Bergmann kernel is concerned) 
to a change of k. Since the W^'^^'s and F(s)'s modular dependence is only in the Bergmann 
kernel, their modular properties are entirely encoded in their k dependence. Thus in [TB] , 
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the following quantities were computed: 



^ and ^ (3.5) 

and they showed a striking similarity with holomorphic anomaly equations. 

The goal of this article is to prove that indeed one can derive the holomorphic anomaly 
equation from the construction of [T6] . 

3.2 Variations of the matrix model free energies 

Let us consider an algebraic curve H{x, y) of genus n and a canonical basis of cycles on it: 

Vi, J = 1 . . . n , AiC] Bj = 5,j , Air]Aj=0 , B^n B^ = 0. (3.6) 

There are n linearly independant holomorphic forms dui on H{x, y) normalized on the A- 
cycles: 

Vz, j = 1 . . . n , (p duj = 6ij (3.7) 
and the Riemann matrix of period r is a symmetric n x n matrix defined by 

dui = Tij. (3.8) 

Bergmann kernel 

There exists a unique bilinear form B_{p, q) with a unique double pole ai p = q without 
residue and normalized on the A cycles: 

+ finite , /b = 0. (3.9) 

where z is any local parameter. B_ is often referred to as the Bergmann kernel in the 
literature. For instance, on a torus, B_ is the Weierstrass function. 

Let us now introduce a new set of cycles, depending on an arbitrary complex symmetric 
matrix n: 

B:=B-tA , A:=A-kB (3.10) 
and define a k— modified Bergmann kernel normalized on these new cycles by the constraints: 

BiP,q) -P-^g , ^^'^!'^^fL + finite , / 5 = 0. (3.11) 
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This definition implies the relation: 

B{p,q) = B{p,q) + 2iiT'^dui{p) Kij duj{q). (3.12) 



Under a modular transformation: 



r 



{C - tD)-^ {tE - F) , E^C-F^D = ld , C^D = D^C , E^F = F^E (3.13) 



the Bergmann kernel 5 changes, as well as the K-modified Bergmann kernel, and the change 
is equivalent to a change of k: 

K D{C - tD)-^ + {{C - TD)-^y k{C - tD)-^ (3.14) 

Therefore, the action of a modular transformation is equivalent to a change of /t, and all 
modular properties of the Bergmann kernel, and thus of the F*^^^'s are encoded in the n 
dependence. 

One can verify that if k = (r — r)~^, then the Bergmann kernel is modular invariant (see 
for instance a proof in [T6]). 

Filling fractions 

We also introduce the filling fractions 

ei := / ydx. (3.15) 



Branch points 

Let us also consider the branch points defined by dx(ai) = 0. When a point p ap- 
proaches a branch point a,, there exists a unique point p such that x{p) = x{p) and p ai. 
Notice that p is defined locally near any branchpoint, but not globally (except for hyperel- 
liptical curves of the form y"^ = P{x)). 

The free energies and correlation functions are recursively defined as follows: 
Correlation functions 

H/^(5)_Q if^<o (3.16) 
Wi^\p) = (3.17) 
wi'\pi,P2) = B{pi,p2) (3.18) 
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and 

Wfc^fi(p,Pi,...,Pfe) 

with the notations (g being near a branchpoint): 

^(g) = -2/®)cia;(g), d£;,(p) = j\{i,p) (3.20) 

where the integration path hes entirely in a vicinity of Oj. If J = {ii,i2, • • • lij} is a set of 
indices, we write pj = {pii,Pi2, ■ ■ ■ ,Pij}, and in the equation above we have K = {1,2, . . . , k} 
and the summation over J is among all subsets of K. 

It is important to notice that by construction, all wjf^'s have vanishing ^-cycle integrals: 



^ W^^^'^ = (3.21) 



which implies: 



Free energies. 

For g > I 



Wj^'^ = K(h Wl:'^ (3.22) 

A Jb 



F(^) = Yl mwi'\q). (3.23) 

^ i 

where $(g) = p ydx is any antiderivative of ydx, i.e. 

d* = ydx (3.24) 

remark: it is proved in that 

^ i 

SO that F^^') can be reg arded as F^^) = ]y(9) _ 
Modular transformations and k dependence. 
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From their definitions, it is clear that the only modular dependence, and also the only k 
dependence of the wj^^^^s is through the Bergmann kernel, and therefore, a modular trans- 
formation is equivalent to a change of n. Moreover, each W^^^ is a polynomial in k, of degree 
at most 3g — 3 + 2k. In order to study modular transformations of the wj^^^^s, or their k 
dependence, it is convenient to compute dW^f^ /dn. This was done in [16] and the result is: 

+ w^-lIm/l 



(3.26) 



and in particular for = 0: 

2ivr-^F(^) = i / / w!f-^\r,s) + \y^l w[''\r) i w'f-''\s) g>2. 

(3.27) 

Also, notice that if k = (r — r)^^, then all wjf^^s and F^^^'s (with 2g + k > 2) are modular 
invariant. 

3.2.1 Variations wrt filling fractions 

The variation of the Bergmann kernel wrt filling fractions follows from Ranch variational 
formula, and it was computed in pTO] that (in fact in [16j, k was considered a constant, and 
thus we just have to add an extra dn/de term): 

dBiv.q) ^ B(r,p)B(r,q)du(r) , ^ dr ^ . . .dn , , . 

)^ = Res ' , J , - 2i7rdu* p /€— Kdu g + 2m du\p)—du{q) 3.28 

oe r-^a dx{r)dy{r) oe oe 

which was written using an operator in |16j : 

--*-(4-|)|: 

SO that: 

B(r,p)B(r, q)du(r) 
A,B{p, q) = Res ^ ^ ' (3.30) 

r-^a dx[r)dy[r) 



This property was sufficient to ensure in |T6], that: 

AMi'^ = -i H^ifi (3.31) 
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4 Holomorphic anomaly equations and matrix models 

In this section we give two proofs of the holomorphic anomaly equation for matrix models. 
In the following, we consider the choice 

1 



T — T 



(4.1) 



in which case the F^^^'s and s are modular invariant [TB] but not holomorphic. 
4.1 A direct proof 

We first consider the variation with respect to e. The only antiholomorphic dependence of 
on e comes from its dependence on k = l/(r — r), therefore 

dW^ _dKdW^ dn__^dT^ ^!^_2i7rr (42) 
de de dn ' de de ' de'^ 

so that 



Then, using eq. (13.261) : 



+ E E ^'tUP'- rWttUPm. ^)) (4.4) 



h LcK 

Notice that with k = l/(r — r), which satisfies: 

= - (4.5) 

the differential operator A of eq. (13.301) reduces to the usual derivative A^^ = d/dej , and 
therefore using eq. (I3.3ip : 



de 



(g) 



Jb Jb J a 



(4.6) 
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Acting again with we find (we use 13.221) : 



and therefore: 



<l + f <\ (4.7) 

+ ir'^^r^ = f f (4-8) 



de^ de de JsJb 
Finally we get the equation: 



de (2i7r)2 de^ 2 V de'^ ' de^ de 



h LCK 

If we normalize the moduli as: 



= (2i7r)^e^ (4.10) 



that last equation reads: 



dt of 2\ dt^ dt dt 

+ E E -d^ — h^) (4-11) 



dt dt 

h LcK 



For the free energy, we set = 0: 



+ ^ (4-12) 



at af 2 V 'at dt (9t at 

Since 

-/t^^ _ /s:-^^ = C|/, ^«:^^^ = -rf, (4.13) 

at^^at^at^ ^' at^ 

we recover the holomorphic anomaly equation fl2.9p for the free energies. The holomorphic 
anomaly equation for the correlation functions can be written (with the convention that 
DiWj^^ = dui) as 

^Wlf^ = (D:DjWt'' + E E DjWi'^DjWt,'^) . (4.14) 

h LcK 
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In the context of topological string theory, these correlation functions correspond to open 
string amplitudes (see [25] for explicit results in the case of mirrors of toric CY manifolds), 
therefore the above equations can be regarded as holomorphic anomaly equations for open 
string amplitudes. These equations have been discussed in a general context in P, H] , and 
they are obtained by looking at the degenerations of Riemann surfaces with boundaries. The 
result (14.141) seems to indicates that the only degenerations that contribute in this case are 
the ones which involve closed string intermediate states. The first term in the r.h.s. of (14.141) 
corresponds to a node, while the sum in the second term corresponds to all the possible 
splittings of the Riemann surface Hg^k into two Riemann surfaces T^hj and Ilg-h,k-i- 
It is also easy to show that the correlation functions (12.121) are given by 

Ct-i. = 7^i <'\pk) (4.15) 

(2i7r)2 Jbi^ Jb,^ 

since the covariant derivative Dj = di — Tj commutes with when acting on a differential 
f{p)dp with vanishing ^-cycle integrals: 

Dj^J{p)dp = j^Dif{p)dp + Di(^j^f{p)dp 

= f Dif{p)dp -Ti (b f{p)dp - diT (b f{p)dp 

JB JB J A 

= <h Djf{p)dp-Ti (h f{p)dp - diT K f f{p)dp 

JB JB JB 



j Dif{p)dp 



(4.16) 



Note that the last equality requires explicitly the condition 

i f{p)dp = i f{p)dp -kI f{p)dp = 0. (4.17) 
J A J A Jb 

If we now integrate (I4.14p over the cycles Bj^ Bj^, , one finds the holomorphic anomaly 
equation for the correlation function (I2.13P in the local case (i.e. without the last line). This 
indicates again that the open string anomaly equations (I4.14p involve only the closed string 
sector in the intermediate states. 



4.2 A combinatorial proof 

In the previous section we have given a simple proof of the holomorphic anomaly equations. 
In pp (see also [I^) it was noticed that these equations have a canonical solution of the form 

F(^)(t,t) = F(^)(t) + r,((Imr)-\9,, ■ ■ ■9,„F,<,), (4.18) 
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where Tg is a polynomial functional in (Imr)~^ and in the derivatives of lower genera holo- 
morphic free energies. In [1], the coefficients of Tg were interpreted diagrammatically as the 
possible degeneracies of a genus g Riemann surface. In this canonical solution, the anti- 
holomorphic dependence of F^^\t,t) enters only through Imr. We will now show that this 
canonical solution can be obtained by iterating f l3.26p . This can be used to obtain a second 
proof of the holomorphic anomaly equations. 



k I k-l 




Figure 1: A graphic representation of the equation fl3.26p . 

Let F^^\t, K = 0) he the holomorphic free energy of the matrix model at genus g. When 
we turn on k we get the K-dependent function 



3.-3 ^ ^ 



F^^\t,K) = J2 '^™^T^^^'ni,0). (4.19) 

ml an™- 



m=Q 



The successive derivatives of F^^\t, k) at k = can be easily computed by iterating (13.26^ . 
We take into account that at k = we can substitute the contour integrals of the correlation 
functions by ordinary {not covariant) derivatives of F^^\ as it follows from fl4.7p . It is 
convenient to represent fl3.26p graphically, so that the terms obtained in the iteration are 
represented by graphs. This representation is shown in Fig. [1] Since the number of legs 
generated in this way is always even, we have labeled the connecting legs by 2z — 1 and 2i. 
The result of the iteration is simply (this is easily seen from the graphical representation in 
Fig. DP): 

Sg— 3 n 

E^G., (4.20) 

m=0 /l,...,/277i = l ' Gm 

where Gm is a connected degenerate surface with m propagators connecting 2m points la- 
belled by 1, . . . , 2m, and in such a way that the point labeled by 2i — 1 is connected by a 
propagator to the point labeled by 2i, z = 1, ■ ■ ■ , m. Each of these surfaces leads to a term 
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Ag^- These terms are constructed as follows: if Gm is made out of r Riemann surfaces of 
genera gi, - ■ ■ ,gr with ni, ■ ■ ■ , punctures, respectively, (such that X] '^j = 2m), and if the 
rii punctures of the i^^ Riemann surface are labeled Ji = { . . . , Ji,ni}, then 



(4.21) 



i=l 



As an example of this procedure, we show in Fig. [2] the graphs that contribute to -^2(^5 1^)- 



i^2W = F2(0) + 




+ 



22-2! 



'1 '2 

^3, 4i . 



23-3! 



+ 
+ 
+ 
+ 








—6 









+ 














— 6 




Figure 2: The graphs that contribute to F2{t, k) after iterating Fig. [H 
Combinatorics 
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Notice because of the sum over indices Ii,...,l2m in (I4.20p . many graphs Gm's give 
the same contribution to F^^\t,K), and thus they better be computed only once with a 
multiphcity factor. The muhiphcity factor is thus A/G'^/(2™m!), where A/g^ is the number 
of inequivalent ways of relabehng the 2m punctures so that the puncture 2? — 1 is hnked by 
a propagator to 2i for alH = 1, . . . , m. Then notice that: 

' (4,22) 



2'^m! 2m! 

and notice that (2m — 1)!!A/g!„ is the number of inequivalent ways of relabeling the 2m 
punctures (without any constraint on the labels of punctures which are linked together), 
i.e. it is the number of inequivalent pairings (compatible with Gm) of 2m points. This 
is exactly the multiplicity factor arising from Wick's theorem, i.e. F^^\t,K) is exactly the 
Feynman graph expansion of an integral whose propagator is N~'^k and whose vertices are 
the iV2-25« Qj^^ ^■■■djj^^ F(3-)(t, 0). 

Notice that because we have chosen W^^^ = and W^i"'' = 0, A.^ vanishes if Gm 
contains surfaces of genus zero with less than 3 punctures. Therefore there is no vertices 
corresponding to the first and second derivative of F^^\ 

Notice also that the recursive application of Fig. [T] generates only connected graphs Gm's, 
and thus we are in fact computing the Feynman graph expansion of the log of an integral. 

Therefore 

oo 

F{t, k) = J2 N^'^'F'^s^t, k) = In {Z{t, k)) (4.23) 

9=0 

is exactly the Feynman graph expansion the following integral 

Z{t,K)= [ dr]exp\--N\ri - t)K'\r] - t) + F{ri,0) - N^ri - tydiF^^\t,0) 

^ 2 (4.24) 

--N\v-tydijF('\t,o){v-ty}. 

Indeed, the Feynman graph expansion, is nothing but the saddle point approximation, and 
here the saddle point equation to large leading order is Kjj{r] — t)"^ — diF^'^\r],0) + 
diF^^\t, 0) + d]jF^'^\t, 0){r] - tY = 0, which is clearly satisfied atr] = t. The propagator is 
simply 

and the vertices are given by 

Ar2-29 

-^dj,---dj,F^^\t), (4.26) 
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where k > 3 ioi g = and k > 1 for g > I. 

We will now show that, for the choice k, = (r — r)^^, the generating functional (14.241) 
satisfies the holomorphic anomaly equations of [6j in the local case. We first note that there 
is a cancellation, 

Z{t,{T-r)-') = j dr/exp|-^iV2(r/-t)f(r/-t)+F(?7,0)-Ar2(r^-t)^9/F(°)(t,0)} (4.27) 
We compute now: 

^diZ = {fjj - Tjj) {{V - tY) + N'diF'~'\ (4.28) 
where ( ) denotes an average in the integral (14.271) . It follows that, 

^djZ = N' {fu-Tu){{v-ty), (4.29) 

Zj 

and 

= -iV2(f,,-r,,)-iV2ajr,x((r7-t)^)+iV^(r-f),K(r-f),,.((r^-t)^(r^-t)^). (4.30) 

Zj 

On the other hand, we compute: 

= Af''CTj^^{^ - tYiv - tV) + Idj^A (4.31) 

where 

CjjK = ^iJ- (4.32) 

and 

A = log det (4.33) 
and we used that the non-holomorphic F^^\k) generated by the integral (I4.27P is 

F^^\t,K) = F^^\t) + A{K). (4.34) 

We express now the quadratic average in (I4.3ip in term of the double derivative (I4.30p . If 
we use 

= K^'dj^fij, (4.35) 

we finally find that 

N"^ ~ 1 / - -\ 

-jdj^Z = -^C'-i [d]jZ - TfjdLZ) , (4.36) 
20 



where the tensor C-^ and the Christoffel symbol Tfj are given by fl2.16p and (12.171) . re- 
spectively. In other words we obtain the holomorphic anomaly equations (12.211) which is 
equivalent to (I2.9p . 

The expansion in [T] is precisely the one that is obtained by expanding the integral (I4.24p . 
This means that the iteration of Fig. [1] leads to a representation of F^^\t, k) which agrees 
with the canonical solution of [Tj for the choice ( 14. ip . 
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